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Abstract: 


[Objective] To study the main differential metabolites of Xihuang Pill (XHP) in rat serum and the mechanism of 
related pathways of metabolites on breast cancer. 

[Methods] Qualitative and quantitative analysis of metabolites in XHP drug serum group and blank serum group 
by liquid chromatography-mass spectrometry (LC-MS) technology. And using sample correlation heat map and 
multivariate statistical analysis methods to compare the metabolic differences between the two groups. The 
metabolites were analyzed by cluster analysis, Variable Importance in Projection (VIP) analysis, Kyoto 
Encyclopedia of genes and Genomes (KEGG) compound classification and KEGG functional and enrichment 
topology analysis. 

[Results] LC-MS technology identified a total of 765 metabolites in the XHP drug serum group; a total of 697 
metabolites in the blank serum group. VIP analysis screened the top 30 serum differential metabolites that were 
significantly different between the two groups, such as Abscisic acid, Quillaic acid, 
2,2-Bis(4-hydroxyphenyl)-1-propanol, Corey PG-Lactone Diol, (S)-Naproxen and so on. KEGG compound 
classification showed that most of the metabolites in XHP were classified as phospholipids and amino acids, 
steroid hormones and carboxylic acids . KEGG functional pathways main involved are Lipid metabolism, Amino 
acid metabolism, Cancer: overview. KEGG enrichment and topology analysis, mainly involved in steroid hormone 
biosynthesis pathway and beta-alanine metabolism pathway. 

[Conclusions] The main differential metabolite of XHP in rat serum may be Abscisic acid. XHP may exert its 
pharmacological effect on breast cancer by regulating steroid hormone biosynthesis pathway to regulate estrogen 
and progesterone levels and beta-alanine metabolism pathway to induce cancer cell apoptosis. 
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1. Introduction 

Xihuang Pill (XHP) is from the “waike zhengzhi quansheng ji" compiled by Wang Weide in the 
Qing Dynasty. Its prescription consists of NIU HUANG, SHE XIANG, RU XIANG, and MO 
УАОГ!. As a classical prescription, XHP has been used in the clinical treatment of various diseases. 
XHP has good curative effects on breast cancer, liver cancer, gastric cancer, colon cancer, lung 
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cancer and other tumors). Modern pharmacological research shows that XHP has antitumor, 
anti-inflammatory, antibacterial, and immune-enhancing effects?! Experiments have demonstrated 
that XHP can restrain the proliferation of tumor cells while promoting tumor cell apoptosis, 
preventing tumor invasion and metastasis, preventing tumor angiogenesis, and regulating the 
tumor microenvironment!^?!, In the early stage, we mainly studied the mechanism of XHP in 
treating precancerous lesions of breast cancer. We used 7,12-dimethylbenzanthracene (DMBA) 
combined with estrogen and progesterone to induce a rat model of precancerous lesions of breast 
cancer and intervened with XHP. XHP may advance the apoptosis of precancerous cells and 
decrease the hyperplasia of vascular endothelial cells by regulating the expression of genes and 
proteins that have a bearing on the PIBK/AKT/mTOR pathway, thereby inhibiting the progression 
of precancerous lesions of breast cancerl°l. XHP may treat precancerous lesions of breast cancer 
by improving microcirculation in rats with precancerous lesions of breast cancer"!. XHP can block 
and reverse the histopathology of breast tissue in rats with breast cancer precancerous lesions 
sequentially triggered by DMBA combined with estrogen and progesterone, and its mechanism 
may be related to its regulation of the expression of NF-kB protein in breast tissuel?!. In clinical 
studies, XHP combined with chemotherapy to treat breast cancer can enhance clinical efficacy, 
improve the quality of life, and reduce the untoward reactions caused by chemotherapy"! 

Metabolomics is used to detect the dynamic and overall changes in metabolites, obtain meaningful 
biomolecular information, and understand the biological process of the body under the actions of 
drugs!!°l, Metabolomics helps to understand the relationship between traditional Chinese medicine 
(TCM) and diseases from the perspective of molecular biology and to study their mechanisms of 
action, which is conducive to the better efficacy of TCM in the body. Metabolomics is now widely 
used in the screening of active components of TCM and the study of the pharmacokinetics and 


1]. Due to 


mechanisms of action of TCM, providing a new method for the modernization of TCM! 
the multicomponent, multitarget and multimechanistic treatment of diseases by the Chinese 
medicine XHP, its complexity and the diversity of processes in the body make it difficult to 
deepen modern research on XHP, and the metabolites of XHP are not yet clear. To solve this 
problem, this article mainly applies LC-MS technology to carry out nontargeted metabolomic 
analysis of the components of XHP. To explore the main differential metabolites of XHP, 
metabolites were classified, and their relationships with diseases were evaluated. Studying the 
differential metabolites of XHP is of great significance for a comprehensive understanding of their 
pharmacological mechanisms in the treatment of breast cancer. 

2. Materials 

2.1 Instruments and Equipment 

NewClassic MS electronic balance (NewClassic MF MS105DU, METTLER TOLEDO); Freezing 
centrifuge (Centrifuge 5424 R, Eppendorf); Temperature Control Ultrasonic Cleaner-10L 
(SBL-10TD, Ningbo Xinzhi Biological Technology Co., Ltd.); Multi sample freeze grinder 
(Wonbio-96c, Shanghai Wanbai Biotechnology Co., Ltd); Bench type rapid centrifugal 
concentration dryer (LNG-T88, Taicang Huamei biochemical instrument factory); Nitrogen 
purging instrument (JXDC-20, Shanghai Jingxin Industrial Development Co., Ltd); UHPLC liquid 
chromatography system (Vanquish Horizon system, Thermo Scientific); Mass spectrometer 
(Q-Exactive, Thermo Scientific). 

2.2 Drugs and Reagents 

XHP (Zhejiang Tianyitang Pharmaceutical Co., Ltd. Division, lot No.: 1703011); Methanol 


(Fisher Chemical); Acetonitrile (Fisher Chemical); Formic Acid (CNW); Water (Fisher Chemical); 
2-Propanol (Merck); 2-Chloro-L-Phenylalanine (Adamas-beta). 

2.3 Animals and Groups 

Twelve SPF 6-week-old female Sprague-Dawley (SD) rats (weighing 180 + 20 g) were purchased 
from Hebei Experimental Animal Center (Shijiazhuang, Hebei, China, license number: 1705351). 
Animal experiments were approved by the Medical Ethics Committee of Hebei University of 
Chinese Medicine and were carried out in accordance with the Animal Welfare Guidelines of the 
Medical Ethics Committee of Hebei University of Chinese Medicine (ethical review number: 
DWLL2019036). The rats were divided into two groups: an XHP drug serum group and a blank 
serum group, with six rats in each group. The rats in the two groups were allowed to feed and 
drink freely, and the experiment was carried out after 1 week of routine feeding. 

3. Methods 

3.1 Preparation of Serum 

3.1.1 Preparation of XHP Drug Serum 

The rats in the XHP-containing serum group were given XHP Chinese medicine solution every 
day. XHP was ground, crushed and made into a suspension with distilled water. The concentration 
was 0.25 g/ml, and the administration volume was 5 ml/kg. It was administered twice a day. The 
cumulative daily dose was 2.5 g/kg. Dose conversion between humans and rats was carried out 
with reference to the “Methodology of Pharmacological Research of Traditional Chinese Medicine” 
edited by Chen Qi. It was administered by gavage according to the dose for 3 days. One hour after 
the last administration on the 4th day (after fasting for 12 h prior to gavage), the rats were 
anesthetized with pentobarbital sodium 200 mg:kg", ip, blood was collected from the abdominal 
aorta, and serum was separated in a sterile tube. After the serum was inactivated at 56 °C for 30 
min, it was stored at -80 °C for later use. 

3.1.2 Preparation of Blank Serum 

The rats in the blank serum group were administered distilled water ig every day (the dosage 
volume was the same as that in the XHP drug serum group). For 3 consecutive days, blood was 
collected and separated at the same time as that in the XHP drug serum group. The method of 
serum preservation was the same as that in the XHP drug serum group. 

3.2 Serum Grouping 

Two hundred microliters of serum from each rat in the XHP serum and blank serum groups was 
extracted for LC-MS metabolomics analysis. The drug serum extracted from the XHP-containing 
serum group was set as the experimental group (sy), and the sy group was divided into 
experimental Group 1 (syl), experimental Group 2 (sy2), experimental Group 3 (sy3), 
experimental Group 4 (sy4), experimental Group 5 (sy5) and experimental Group 6 (sy6). The 
serum of the blank serum group was set as the control check (ck), and ck group was divided into 
control check 1 (ck1), control check 2 (ck2), control check 3 (ck3), control check 4 (ck4), control 
check 5 (ck5) and control check 6 (ck6). The above experimental group included 6 groups, and the 
control group included 6 groups. LC-MS nontargeted metabolomics analysis and multivariate 
statistical analysis were performed to compare the metabolic differences between the two groups. 
3.3 LC-MS Analysis 

3.3.1 Sample Treatment 

We accurately removed 100 uL serum from each group to a 1.5 mL centrifuge tube and added 400 
uL extract (methanol:acetonitrile = 1:1 (v:v)) containing 0.02 mg/mL internal standard 


(L-2-chlorophenylalanine); after vortex mixing for 30 s, we performed low-temperature ultrasonic 
extraction for 30 min (5 °C, 40 kHz), placed the sample at -20 °C for 30 min, centrifuged for 15 
min (13000 x g, 4 °C), took supernatant and dried the sample with nitrogen, added 120 uL 
reconstituted solution (acetonitrile:water = 1:1) to reconstitute, vortexed and mixed for 30 s, 
followed by low-temperature ultrasonic extraction for 5 min (5 °C, 40 kHz), centrifuged for 10 
min (13000 x g, 4 °C), and transferred the supernatant to an injection vial with an inner cannula 
for analysis. In addition, 20 uL of the supernatant was pipetted for each sample and mixed as a 
quality control sample. 

3.3.2 Chromatographic Conditions 

The chromatographic column was an ACQUITY UPLC HSS T3 (100 mm x 2.1 mm i.d., 1.8 um; 
Waters, Milford, USA); mobile phase A was 95% water + 5% acetonitrile (containing 0.1% formic 
acid), and mobile phase B was 47.5% acetonitrile + 47.5% isopropanol + 5% water (containing 
0.1% formic acid). The flow rate was 0.40 mL/min, the injection volume was 2 uL, and the 
column temperature was 40 °С. 

3.3.3 Mass Spectrometry Conditions 

The sample was electrospray ionized, and mass spectrometry signals were collected in positive 
and negative ion scanning modes. The scan type parameter was 70-1050 m/z; the sheath gas flow 
rate parameter was 40 arb; aux gas flow rate parameter was 10 arb; the heater temp parameter was 
400 °C; the capillary temp parameter was 320 °C; the spray voltage (+) parameter was 3500 V; the 
spray voltage (-) parameter was -2800 V; the s-Lens RF Level parameter was 50; the normalized 
collision energy was 20, 40, and 60 eV; the resolution (Full MS) was 70000; and the resolution 
(MS2) was 17500. 

3.3.4 Preparation of Quality Control Samples 

Quality control (QC) samples were prepared by mixing the extracts of all samples in equal 
volumes. The volume of each QC was the same as that of the sample. It was processed and tested 
in the same way as the analysis sample. QC was injected every 4 samples to check the stability of 
the whole detection process. Therefore, the quality control samples of this experiment are divided 
into quality control sample 1, quality control sample 2, and quality control sample 3. 

3.4 Data Analysis 

Thermo Fisher's Ultra High-Performance Liquid Chromatography tandem Fourier Transform mass 
spectrometry UHPLC-Q Exactive system was used for chromatographic peak identification. 
Progenesis QI software (Waters Corporation, Milford, USA) was used to perform baseline 
filtering, peak identification, integration, retention time correction, and peak alignment to obtain 
information such as the mass-nucleus ratio (m/z), retention time and peak area. Then, the software 
was used to search and identify the characteristic peaks, and the MS and MS/MS mass spectral 
information was matched with the metabolic database. The MS mass error was set to be less than 
10 ppm, and the metabolites were identified according to the secondary MS spectrometry 
matching score. The main databases were http://www.hmdb.ca/, https://metlin.scripps.edu/and 
additional mainstream public databases and self-built databases. The normalized data matrix was 
imported into the Ropls software package for pattern recognition analysis of PCA and OPLS-DA. 
Analysis software platforms and related versions involved in this paper see Table 1. 

4. Results 

4.1 Venn Analysis of Differential Metabolites 

The substances detected by the mass spectrometer were obtained through LC-MS analysis 


(Supplementary Table S1). Through primary and secondary mass spectrometry data, libraries 
(self-built database, Metlin database, HMDB database, etc.) were searched to finally obtain the 
annotated substances. See Supplementary Tables S2 and S3 for the metabolite information of the 
experimental and control groups, respectively. A total of 765 XHP metabolites were identified in 
the serum of the experimental group; 697 serum metabolites were identified in the control group. 
The obtained metabolite data of the experimental and control groups were imported into Venn 
Diagram (R package) software, and the intersection was taken to obtain the coacting metabolites. 
Venn diagram analysis revealed the number of metabolites and the overlapping relationships of 
metabolites in each group. The results showed that among the positive ion metabolites, there were 
396 in the control group and 442 in the experimental group, and the intersection of the two 
positive ion metabolites yielded 396; among the negative ion metabolites, there were 301 in the 
control group and 323 in the experimental group, and the intersection of the two negative ion 
metabolites yielded 301. The results are summarized in Figure 1. 

4.2 Comparison and Screening of Two Groups of Differential Metabolites by 

4.2.1 Difference Statistics and Volcano Map 

The two groups of differential metabolites were compared and graphed, as shown in Figure 2. The 
X-axis is the multiple shift value of the metabolite expression difference between the two groups, 
1.е., log2FC, and the Y-axis is the statistical test value of metabolite expression difference, 1.е., 
-10210 (p value). The higher the value is, the more remarkable the expression discrepancy is. The 
values of the X- and Y-axes are logarithmic. Every point in the figure represents a given 
metabolite, and the point size represents the VIP value size. The point on the left is the metabolite 
of differentially downregulated expression, and the point on the right is the metabolite of 
differentially upregulated expression. The closer it is to the left, right and upper points, the more 
remarkable the expression discrepancy. The results showed that the experimental and control 
groups were significantly different. There were 115 positive ions and 87 negative ions of 
differential metabolites, and the expression levels of differential metabolites were upregulated 
more than downregulated. See A, B and C in Figure 2. 

4.2.2 Details of the Differences Between Groups 

The differential metabolism between the groups was studied, and a histogram of the relative 
expression abundance of metabolites in each group of samples was obtained. The difference test 
methods included the Student's t test (Unpaired) using a two-tailed test; OPLS-DA data conversion 
was performed using Pareto transformation, for which the OPLS-DA confidence level was 0.95, 
and the OPLS-DA number of permutations was 200 (P value < 0.05, VIP pred OPLS-DA > 1, 
display FC value: > 1 or < 1). The Y coordinate is the mass spectrum intensity value (the mass 
spectrum intensity after data preprocessing), and the error bars represent the average 
value+standard deviation. The VIP value was selected according to OPLS-DA in descending order, 
and the largest metabolite was selected for differential metabolism analysis, as shown in D, E and 
F in Figure 2. Abscisic acid is the metabolite with the largest difference between the experimental 
and control groups in positive ion mode and positive-negative ion mode. In negative ion mode, 
(S)-naproxen was the metabolite with the greatest difference between the experimental and control 
groups. The contents of these two metabolites in the experimental group were greater than those in 
the control group. 

4.2.3 PCA Analysis 

The PCA method was used to analyze the LC/MS data of the experimental and control groups. 


Then, the samples were analyzed by dimensionality reduction. The PCA results in positive and 
negative ion modes are shown in Figure 3, A and B, respectively. The experimental and control 
groups showed obvious separation trends, indicating that the metabolites of the two groups were 
significantly different. The two sets of samples could be separated by PCA. The confidence ellipse 
implies that the "true" samples of this group are distributed in this area with a 95% confidence 
level; more than this area can be considered a probable abnormal sample. Panels A and B in Figure 
3 show that no abnormal samples that deviated from the confidence circle were detected in the 
LC-MS. Green dots, control group; blue triangles, experimental group. 

4.2.4 PLS-DA Analysis 

The LC/MS data of the two groups of samples were analyzed by the PLS-DA method, and the 
samples were analyzed by dimensionality reduction. The PLS-DA analysis results in positive and 
negative modes are shown in Figure 4, A and B, respectively. The experimental and control groups 
showed obvious separation trends, indicating that the metabolites of the two groups were 
significantly different. 

The serum metabolic profiles of the two groups showed significant separation in positive and 
negative ion modes. To test the reliability of the model and prevent overfitting, the model 
overview and displacement test chart were used for verification. The PLS-DA model is shown in 
panels C and D in Figure 4. There were two principal components in the positive ion model: the 
model interpretation rate R2X was 0.407, the predictive ability parameter О? was 0.953, and Q?> 
0.5, indicating that the model had a better predictive ability and a better fit. There were two 
principal components in the negative ion model: the model interpretation rate R?X was 0.523, the 
predictive ability parameter Q2 was 0.966, and Q2> 0.5, indicating an excellent degree of fit. 

In the displacement test, the abscissa represents the displacement retention, the ordinate represents 
the values of the R? (red dots) and Q? (blue triangles) displacement tests, and the two dotted lines 
represent the regression lines of R2 and Q2, respectively. Permutation tests were performed for the 
experimental and control groups. The replacement test model randomly disrupts the grouping 
labels (Y variables) of the experimental and control groups. The replacement retention of the 
horizontal axis was the proportion consistent with the order of the Y variables of the original 
model. The replacement retention was 1, that is, R? and О? of the original PLS-DA model. The 
number of random permutation tests was 200, R? was 0.8818, and Q? was -0.0401 (negative ion 
model). В? was 0.9435, and Q? was -0.0710 (positive ion model). The negative value of Q? 
indicated that the model built was reliable, there was no overfitting phenomenon, and the 
modeling was successful. As the replacement retention decreased, the regression line showed an 
upward trend, indicating that the replacement test was passed (Figure 4, E and F). 

4.2.5 OPLS-DA Analysis 

The OPLS-DA score map filters out information irrelevant to the grouping through orthogonal 
rotation, thereby better making a distinction between groups and improving the function of the 
model. The OPLS-DA results in positive and negative modes are shown in Figure 5, A and B, 
respectively. The experimental and control groups showed obvious separation trends, indicating 
that the metabolites of the two groups were significantly different. 

To verify the reliability of the model and prevent overfitting, the model overview and permutation 
test chart were used for verification. The OPLS-DA model, is shown in panels C and D in Figure 5. 
There were two principal components in the positive ion model: the model interpretation rate R2X 
was 0.407, the predictive ability parameter Q2 was 0.919, and Q? > 0.5, indicating that the model 


had a better predictive ability and a better fit. There were two principal components in the negative 
ion model: the model interpretation rate R2X was 0.523, the predictive ability parameter О? was 
0.96, and Q2 > 0.5, indicating an excellent degree of fit. 

In the displacement test, the abscissa represents the displacement retention, the ordinate represents 
the values of the R? (red dots) and Q? (blue triangles) displacement tests, and the two dotted lines 
represent the regression lines of R2 and Q2, respectively. The number of random permutation tests 
was 200, В? was 0.8811, and О? was -0.0948 (negative ion model). В? was 0.943, and О? was 
-0.0607 (positive ion model). The negative value of Q2 indicated that the model built was reliable, 
there was no overfitting phenomenon, and the modeling was successful. As the replacement 
retention decreased, the regression line showed an upward trend, indicating that the replacement 
test was passed (Figure 5, E and F). 

4.3 Cluster Analysis of Metabolites 

Hierarchical cluster analysis was carried out on the different metabolites in the experimental and 
control groups. The metabolite clustering analysis algorithm was based on hierarchical clustering. 
The subclusters were used select the top ten, and the top 30 most abundant metabolites were 
selected, as shown in Figure 6. On the left is a dendrogram of metabolite clustering, and the 
metabolite name is on the right. The closer the two metabolite branches are, the more similar their 
expression levels are; the dendrogram of sample clustering is at the top, and the name of the 
sample is at the bottom. The closer the two sample branches are, the closer the expression pattern 
of all metabolites of the two samples is, that is, the closer the shift trend of metabolite expression 
is. The metabolite expression trends of the control group and the experimental group were 
different. In the control group, the expression levels of 2-indolecarboxylic acid and 
indole-3-carboxilic acid-O-sulfate in subcluster 5 and carboxyoprost bromethamine in subcluster 8 
and (+)-propionylcarnitine in subcluster 10 were upregulated, while other subclusters were 
downregulated. The expression levels of the remaining metabolites were upregulated in the 
experimental group. 

4.4 Variable Importance in Projection (VIP) Value Analysis 

According to the weighting coefficient of the OPLS-DA model, the VIP score was used to rank the 
identification contributions of the metabolites in the experimental and control groups. The 
metabolite clustering algorithm was hierarchical clustering, the metabolite distance algorithm was 
Euclidean, the metabolite hierarchical clustering method was complete, the VIP value was derived 
from OPLS-DA, VIP > 1, and the metabolites with the top 30 VIP values were selected. As shown 
in Figure 7, a higher VIP value indicated that the metabolite differed greatly between the two 
groups. The top five metabolites were abscisic acid (VIP = 6.1604), quillaic acid (VIP = 6.0746), 
2,2-bis(4-hydroxyphenyl)-1-propanol (VIP = 5.543), Corey PG-lactone diol (VIP = 5.4712) and 
(S)-naproxen (VIP = 4.7948). The selected differential metabolites are shown in Table 2. 

4.5 KEGG compound classification 

Kyoto Encyclopedia of Genes and Genomes (KEGG) compound classification is based on 
compounds with biological roles. The identified metabolites were compared to the KEGG 
Compound database to obtain the metabolite classification profile and statistical mapping. KEGG 
compound classification revealed that the metabolite classification in XHP consisted of mostly 
phospholipids (22 metabolites), amino acids (14 metabolites), steroid hormones (6 metabolites) 
and carboxylic acids (5 metabolites), as shown in Figure 8. 

4.6 KEGG functional pathway analysis 


By comparing the results with the KEGG database, the metabolites in the experimental group can 
be sorted in light of the pathways touched upon or the functions they perform. As shown in Figure 
9, the main pathways involved included lipid metabolism (64 metabolites), amino acid metabolism 
(44 metabolites), cancer: overview (41 metabolites), digestive system (28 metabolites), 
metabolism of other amino acids (16 metabolites), nervous system (15 metabolites), and 
membrane transport (15 metabolites). XHP had the characteristics of multicomponent, multitarget 
and multilink overall effects. To further study the mechanism of XHP in treating breast cancer, we 
sorted out the metabolites (15 numbers) that XHP significantly regulated that were related to 
breast cancer and their functional predictions (see Supplementary Table S4). 

4.7 KEGG Pathway Enrichment Analysis 

Through KEGG enrichment and topological analyses, 13 different metabolic pathways were found, 
as shown in Figures 10 and 11. KEGG enrichment analysis showed that the two groups had 
significant changes in the following pathways: necroptosis, steroid hormone biosynthesis, 
beta-alanine metabolism, vitamin digestion and absorption, and apoptosis. Topological analysis 
revealed 4 main metabolic pathways, namely, steroid hormone biosynthesis, beta-alanine 
metabolism, retinol metabolism, and caffeine metabolism. The metabolic pathways with 
statistically significant differences (P « 0.05) were the steroid hormone biosynthesis pathway and 
the beta-alanine metabolism pathway. The steroid hormone biosynthesis pathway included three 
metabolites: 21-hydroxypregnenolone (c05485), corticosterone (c05488) and 21-deoxycytortisol 
(c05497). The beta-alanine metabolism pathway included two metabolites: carnosine (000386) and 
pantothenic acid (c00864). See Figure 12. 

5. Discussion 

With the gradual deepening of biotechnology research, some molecular compounds have been 
found to have impacts on some diseases. For example, F. S. Mohammed et al.!?! showed that 
sechium edule extract results in the prolongation of QT or QTc and RR intervals and increases the 
heart rate in high-fat diet-challenged mice. S. Taifa et а1.П31 found that copper nanoparticles have 
an excellent effect on mastitis. This article studies the metabolic compounds of XHP and the 
mechanisms of its metabolites in the treatment of breast cancer XHP is composed of NIU 
HUANG, SHE XIANG, RU XIANG, and MO YAO. It is a commonly used Chinese patent 
medicine for the treatment of tumors, especially breast cancer and other diseases. In recent years, 
there have been many pharmacological studies on the four traditional Chinese medicines in XHP. 
These studies have shown that NIU HUANG and SHE XIANG inhibit the proliferation of cancer 


1^ 131 In addition, the active ingredients in NIU HUANG can reduce the proliferation of 


cells! 
human colon adenocarcinoma cells! Deoxycholic acid, a component of NIU HUANG, can 
block cells in G2/M phase, inhibit the growth of human gastric cancer BGC-823 cells and induce 
apoptosis through the mitochondrial pathway!'7!. Muskone, the active ingredient in SHE XIANG, 
can significantly inhibit the expression of basic fibroblast growth factor and vascular endothelial 
cell growth factor in the nude mouse model of breast cancer MDA231 cell line, thereby inhibiting 
tumor angiogenesis!9). RU XIANG and MO YAO can promote the apoptosis of cancer cells!!*!, 
The main antitumor components of RU XIANG are boswellic acids, which have various effects, 
such as reducing the proliferation of cancer cells, promoting their apoptosis, and intervening to 
attenuate the conduction of tumor-related signaling pathways?”|, 

At present, the pharmacological mechanism of XHP based on nontargeted metabolomics has not 
been found. Therefore, this study used nontargeted metabolomics to analyze the metabolites of 


XHP in rat serum. Mass spectrometry (MS) and nuclear magnetic resonance (NMR) are the two 
most widely used techniques for compound identification. In compound analysis, MS and NMR 
are complementary relationships. MS can provide the atomic formula of an analyte, while NMR 
provides the structural parts into which these atoms are organized. For example, NMR can 
distinguish between isobaric compounds and positional isomers. MS can recognize certain 
functional groups, such as sulfate and nitro groups. These complementary capabilities of MS and 
NMR mean that neither detector can be employed alone for comprehensive compound analysis in 
the fields of natural product discovery, metabolomics, and drug metabolite identification?", As an 
analytical platform, NMR is a powerful tool for studying the structure and dynamics of 
biomolecules under physiological conditions??!, NMR can play an essential role in drug discovery 
because of its ability to provide protein-ligand binding information in solution?! NMR is the 
earliest technique to measure metabolites in biological specimens, and this technique has been 
successfully applied in selecting tumor markers, judging prognosis, pathological typing, and early 
diagnosis4!, The advantages of NMR techniques include selective isotopic detection in complex 
mixtures, determination of unknown key structural parameter metabolites, inherent quantification, 
in situ analysis of pathway dynamics from cells to whole organisms, and nondestructive effects”. 
However, NMR has many shortcomings: the spectral resolution is not high, and it does not have 
the function of chromatographic separation; moreover, the derived spectrum is affected by the 
mutual interference of the peaks of various compounds, so the analytical results are limited, and 
its sensitivity is also lower than that of LC-MSP9l, LC-MS can make up for these two deficiencies 
at the same time and can be used to conduct more accurate qualitative and quantitative analyses of 
target compounds in samples. 

Through nontargeted metabolomic analysis, we identified a total of 765 metabolites in the XHP 
drug serogroup and 697 metabolites in the blank serogroup. The two groups were compared by 
multivariate statistical analysis, and the results showed that there were metabolite differences 
between the two groups. In this study, through the comparison of different metabolites between 
the two groups, it was found that abscisic acid and (S)-naproxen were significantly different 
between the two groups, and the contents of the two metabolites in the XHP drug serum group 
were higher than those in the blank serum group and may be the main functional components of 
XHP in blood metabolism and circulation. VIP analysis found that abscisic acid, quillaic acid, 
2,2-bis(4-hydroxyphenyl)-1-propanol, Corey PG-lactone diol and (S)-naproxen were significantly 
different in the two groups. Therefore, abscisic acid and (S)-naproxen may be the main differential 
metabolites between the two groups. Abscisic acid is a plant hormone that exists not only in plants 
but also in humans and animals. In mice, abscisic acid reduces obesity and interferes with insulin 
resistance. Abscisic acid can inhibit the proliferation of various cancer cells and induce apoptosis, 
induce the growth of mesenchymal stem cells and promote the growth of hematopoietic progenitor 
cells, and play a certain anti-inflammatory role in the immune system. Therefore, abscisic acid 
may be a potential drug for the treatment of various diseases", Domestic and foreign studies 
have found that abscisic acid can resist tumors and inhibit the formation of blood vessels around 
tumors. Recent studies have found that abscisic acid can also be produced by mammalian cells and 
can inhibit the proliferation of cancer cells and regulate the dormancy of cancer cells in the bone 
тагго\ 8-30], H. W. Zhao et а1.281 found that abscisic acid can enhance the activity of Caspase-3 in 
cancer cell lines, increase cancer cell apoptosis, slow down the growth of human tongue cancer 
Tca8113 xenografts, and induce differentiation. D. S. Jang et al.P!! found that abscisic acid 


induced an increase in the content of quinone reductase, which could inhibit the growth of Hepa 
1с1с7 hepatoma. 

(S)-Naproxen is a nonsteroidal anti-inflammatory drug that can reduce the synthesis of 
thromboxane, prostacyclin and prostaglandin from arachidonic acid by inhibiting the 
cyclooxygenase subtype and has anti-inflammatory, analgesic and antipyretic pharmacological 
effects’? It is used in the treatment of chronic arthritis, osteoarthritis, ankylosing spondylitis, 
gout, and other conditions'*!, (S)-Naproxen exerts its anticancer adjuvant effect only by inhibiting 
the activation of cyclooxygenase-2. Its inhibition of cyclooxygenase-1 can induce the synthesis of 
prostaglandin E in the gastric mucosa, resulting in adverse reactions such as gastrointestinal tract 
upset??l, Z. J. Li et al.P?! found that cisplatin compounds containing naproxen structures showed 
special antitumor activity, and the in vitro cytotoxicity of cisplatin coupled with naproxen 
structures on the breast tumor cell lines MCF-7 and MDA-MB-231 was stronger than that of 
cisplatin alone. In conclusion, abscisic acid and (S)-naproxen have inhibitory effects on the 
progression of breast cancer to a certain extent, but the specific mechanism needs to be further 
explored. 

The functional KEGG pathways mainly involve lipid metabolism and amino acid metabolism. 
Modern studies have found that lipid metabolism and amino acid metabolism pathways have 
certain inhibitory effects on breast cancer. The regulation of lipid metabolism, such as lipid uptake, 
synthesis, and hydrolysis, is critical for maintaining cellular homeostasis?*°!, Dysregulation of lipid 
metabolism is one of the most prominent metabolic alterations in cancer. Cancer cells utilize lipid 
metabolism for energy, biofilm composition, and signaling molecules required for proliferation, 
survival, invasion, metastasis, and response to tumor microenvironment influences and cancer 


37]. Cancer cells in the tumor microenvironment constantly change nutrient availability 


therapy! 
during tumor progression, utilizing lipid metabolism to support their rapid proliferation, survival, 
migration, invasion, and metastasis. Lipid metabolism is one of the conditions for tumor invasion 
and metastasis by increasing the decomposition of lipids. The lipids of tumor cells are composed 
of adipocytes, and fatty acid derivatives produced by the decomposition of adipocytes can enhance 
the invasiveness of breast cancer!*!, By studying the relationship between menopause and tumor 
development, Y. Dunneram found that women have a decrease in estrogen after menopause, which 
may lead to changes in fat metabolism and promote breast cancer'?l. If a woman's estrogen levels 
drop significantly after menopause, this hormonal change is likely to lead to a disorder of lipid 
metabolism. Therefore, postmenopausal women are more prone to visceral lipid accumulation 
than premenopausal women, and the increase in fatty acid derivatives produced by its 
decomposition can induce an increased risk of cancer'®!, Estrogen receptor-positive breast cancer 
is more dependent on estrogen and lipid metabolism. Therefore, the lipid metabolism pathway 
may be closely related to the development of and changes associated with breast cancer, regulating 
lipid metabolism, inhibiting the decomposition of fat cells, and helping to inhibit the growth and 
metastasis of breast cancer. 

Modern research has found that amino acid metabolism is involved in tumor growth, amino acids 
are essential for the survival of tumor cells, and tumor cells have an increased demand for 
nutrients to promote proliferation and cancer progression. The more amino acids the body ingests, 
the stronger the tumor's proliferation ability. Recent studies have shown that the amino acid 
metabolism pathway is altered in breast cancer and that amino acid transporters influence tumor 


growth and progression. In breast cancer, glutamine is one of the key nutrients, and glutamine 


metabolism is closely related to amino acid transporters". The proliferation of cancer cells 
depends on the synthesis of intracellular nucleotides and the continuous replication of DNA. 
Glutamine, serine, glycine and aspartic acid are amino acid raw materials necessary for nucleotide 
synthesis, and serine is involved in DNA methylation. Arginine, leucine and glutamine act as 
signaling factors to activate the mammalian target of rapamycin signaling pathway. This pathway 
is closely related to tumor growth regulation, so the mechanism by which this pathway regulates 
tumors can be studied from the aspect of amino acids. The proliferation and maintenance of tumor 
cells are dependent on the supply of amino acids to the intracellular space, which is regulated by 
amino acid transporters. The development and changes of tumors are caused by many pathological 
molecular mechanisms induced by various factors, and the related mechanisms involved in amino 
acid metabolism include oncogene activation, related protein expression, cancer cell proliferation, 
cancer cell proliferation and metastasis, immune suppression, and immune escape!) 

KEGG enrichment and topological analysis showed that XHP was mainly involved in the steroid 
hormone biosynthesis and beta-alanine metabolism pathways. The three XHP metabolites 
involved in the steroid hormone biosynthesis pathway in this study were 21-hydroxypregnenolone, 
21-deoxycortisol, and corticosterone. 21-Hydroxypregnenolone is involved in the synthesis of 
steroid hormones and can prevent fatigue and reduce stress]. 21-Deoxycortisol is a precursor of 
steroid hormones, and corticosterone is a steroid hormone and a key substance in the biosynthetic 
pathway of steroid hormones. Steroid hormones are fat-soluble hormones, and according to their 
pharmacological effects, they are divided into two cholesterol derivatives: adrenal cortex 
hormones and sex hormones. Steroid hormones can enter target cells by diffusion or carrier 
transport, and after entering cells, they bind to receptors to form hormone-receptor complexes and 
play roles in regulating gene expression. The synthesis of steroid hormones usually takes place in 
the adrenal cortex, gonads, brain, placenta, and adipose tissue, and two types of enzymes are 
mainly involved in their synthesis: cytochrome P450 enzymes and  hydroxysteroid 
депудгосепазез “1. Steroid hormones play important roles in regulating hormone-dependent 
tumor development, maintaining fluid balance, regulating immunity, and mediating changes in 
body metabolism and the stress response [44]. The secretion and growth of SHE XIANG are 
related to the biosynthesis of steroid hormones'**!. RU XIANG and MO YAO participated in the 
biosynthesis of steroid hormones in the treatment of breast hyperplasia), Steroid hormones play 
prominent roles in the pathogenesis and progression of breast cancer. The steroid hormones 
estrogen and progesterone may induce breast cancer progression by increasing breast cancer stem 
cells", In clinical applications, hormone replacement therapy (HRT) is used to treat sex 
hormone-dependent tumors!*!, The steroid hormone biosynthesis pathway synthesizes estrogen 
and progesterone, and modern studies have shown that estrogen can promote the proliferation of 
estrogen receptor-positive breast cancer cells. Therefore, by inhibiting this pathway, reducing 
estrogen production may inhibit the proliferation of estrogen—progesterone receptor-positive breast 
cancer cells. Whether XHP affects the occurrence and development of estrogen—progesterone 
receptor-positive breast cancer by regulating the steroid hormone biosynthesis pathway is worthy 
of further experimental investigation. 

The two XHP metabolites involved in the beta-alanine metabolism pathway in this study were 
carnosine and pantothenic acid. Modern research has found that beta-alanine is a unique B-type 
amino acid in the body and does not participate in protein synthesis. Cytosine and uracil are 


metabolized to obtain beta-alaninel?l, and the precursor substance beta-alanine is involved in the 


synthesis of coenzyme a and pantothenic acid and combines with histidine to form carnosine and 


its derivative апѕегіпебо 51. 


Therefore, in the beta-alanine metabolism pathway, carnosine, 
pantothenic acid and beta-alanine are closely related. Carnosine helps to inhibit tumor 
proliferation, and with the accumulation of carnosine, tumor growth is significantly inhibited. The 
reason is believed to be that the accumulation of carnosine hinders glycolysis and energy 
production'?!, Beta-alanine can transmit information between neurons and can also regulate 
hormones. It is an intermediate metabolite of various active substances?! and regulates various 
physiological functions, such as metabolism in the body°!. Beta-alanine metabolism may be 
involved in the apoptosis pathway. B. Saunders et al.P*! found that in the early stage of apoptosis, 
the continuous increase in beta-alanine content can cause oxidative stress damage to the brain, but 
it can increase the level of carnosine in muscles and relieve soreness. T. Gemelli et al.°! showed 
that the gradual increase in the content of beta-alanine in the body will reduce the activities of the 
two main enzymes in glycolytic enzymes, pyruvate kinase and lactose dehydrogenase, thereby 
destroying the glycolysis process and inducing apoptosis. Y. Xu showedP9! that XHP can regulate 
beta-alanine, inhibit the expression of lactate dehydrogenase in breast cancer 4T1 cells, reduce the 
increase in cancer cells, and promote cell apoptosis. Therefore, XHP may induce breast cancer cell 
apoptosis by regulating the beta-alanine metabolism pathway. 


6. Conclusion 

In conclusion, this study used LC-MS technology to study the difference between the metabolites 
in the XHP and blank serum groups. After comparative analysis, a total of 30 different metabolites 
were obtained. Abscisic acid and (S)-naproxen were the main differential metabolites. Most of the 
metabolites were classified into phospholipids and amino acids, steroid hormones and carboxylic 
acids. The major pathways involved in the KEGG functional pathway were lipid metabolism, 
amino acid metabolism, and cancer. Thirteen KEGG enrichment pathways were identified, and the 
steroid hormone biosynthesis pathway and the beta-alanine metabolism pathway were the main 
pathways involved in tumor proliferation, differentiation and apoptosis. The metabolic changes of 
XHP in rats can elucidate the nature of the main differential metabolites and the main metabolic 
pathways involved, allowing a more in-depth study of XHP in the treatment of breast cancer and 


other tumors. 
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Figure 1 Veen diagram of differential metabolites. 


Note: A. Positive ion mode; B. Negative ion mode. 
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Figure 2. Two groups of difference volcano plots and abundance histograms. 


Note: A. Positive ion mode; B. Negative ion mode; C. Positive-negative ion mode; D. Positive ion mode; E. Negative ion mode; 


F. Positive-negative ion mode. 
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Figure 3 Two sets of PCA score chart. 


Note: A. Positive ion mode; B. Negative ion mode. 
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Figure 4 Two groups of PLS-DA score diagrams, model overview diagrams, and permutation test diagrams. 


Note: A. Positive ion mode score diagram; B. Negative ion mode score diagram; C. Positive ion mode model overview diagram; 
D. Negative ion mode model overview diagram; E. Positive ion mode replacement inspection diagram; F. Negative ion mode 


replacement inspection diagram. 
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Figure 5 Two groups of OPLS-DA score diagrams, model overview diagrams and permutation test diagrams. 


Note: A. Positive ion mode score diagram; B. Negative ion mode score diagram; C. Positive ion mode model overview diagram; 
D. Negative ion mode model overview diagram; E. Positive ion mode replacement inspection diagram; F. Negative ion mode 


replacement inspection diagram. 
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Figure 6 Metabolite cluster analysis chart. 


Note: The color in the figure represents the relative expression level of the metabolite in the group of samples. The specific 


expression level change trend is shown in the number label under the color bar at the bottom right. 
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Figure 7 VIP analysis chart. 


Note: On the right is the heat map of metabolite expression, each column represents a sample, and the sample name below; Each 


row represents a metabolite, and the color represents the relative expression level of the metabolite in the group of samples. 
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Figure 8 KEGG compound classification diagram. 
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Figure 9 KEGG functional pathway diagram. 
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Figure 10 KEGG metabolic pathway analysis diagram. 
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Figure 11 KEGG topology analysis path diagram. 
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Figure 12 Network diagram of KEGG pathway in steroid hormone biosynthesis and beta-alanine metabolism. 


Note: Metabolites on a white background indicate metabolites in this pathway; pathways on a red background indicate 
metabolites in this metabolic concentration. 


Tables 


Table 1. Analysis software platforms and related versions involved in this paper. 


Analysis item 

Metlin database 

HMDB database 

Venn Diagram 

Expression data preprocessing 


Correlation analysis 


Differential Metabolite Analysis - Multivariate 


Statistics 
PCA analysis 
Cluster analysis 


Heatmap 


VIP analysis 


KEGG compound classification 


KEGG pathway enrichment 


Software/URL 
https://metlin.scripps.edu/ 
http://www.hmdb.ca/ 

Venn Diagram (R packages) 
Majorbio's own software 


Scipy (Python) 


Ropls (R packages) 


Ropls (R packages) 
Scipy (Python) 


Scipy (Python) 

Ropls (R); 

scipy (Python) 
https://www.kegg.jp/kegg/compound/ 
Scipy (Python) 


Version 

no 

Version 5.0 
Version1.6.20 
no 


Version1.0.0 


Version1.6.2 


Version1.6.2 
Version1.0.0 


Versionl.0.0 
Version1.6.2; 
Versionl.0.0 
kegg v 94.2 

Versionl.0.0 


ID 


pos_1413 
neg 3106 
neg 2250 
pos 2847 


pos 1437 


neg 7754 
pos 1615 
pos 646 

pos 1583 


neg 2304 


neg 7701 
pos 1381 
pos 643 

pos 673 

pos 2977 
pos 1604 
pos 1490 
pos 2481 
pos 2740 
neg 3797 
neg 447 

neg 3897 
neg 7071 
neg 584 

pos 1426 
pos 3000 
pos 3055 
pos 3212 
neg 7818 


neg 2906 


Table 2 Serum Differential Metabolism Markers of XHP 


Metabolite 


Abscisic acid 
11,12,13-Trinor-1(10)-spirovetivene-2,7-dione 
8-Acetylegelolide 

21-Deoxycortisol 


2,2-Bis(4-hydroxyphenyl)-1-propanol 


13-[2-(3.4-dihydroxy-5-methoxyphenyl)ethyl ]phenylj oxidanesulfonic acid 
Quillaic acid 

Corey PG-Lactone Diol 

PGK2 


1[3-(2-hydroxyphenyl)oxiran-2-yl]methoxy sulfonic acid 


Citreovirenone 

Helenalin 

Glandulone B 

AUSTRICINE 

3-(1,1-Dimethylallyl)scopoletin 
lalpha-hydroxy-22-0xo-23,24,25,26,277-pentanorvitamin D3 
Pterosin E 

24-Norursa-3,12-dien-11-one 

Isoglabrolide 

Glyyunnansapogenin B 

Formononetin 7-(6"-malonylglucoside) 

Fusidic Acid 

Lactapiperanol D 

8-Deoxy-11,13-dihydroxygrosheimin 

Achimilic acid 
(+)-cis-5,6-Dihydro-5-hydroxy-4-methoxy-6-(2-phenylethyl)-2H-pyran-2-one 
2-Butanone, 4-(6-hydroxy-2-naphthalenyl)- 

(4S, 6R)-p-Mentha-1,8-diene-6,7-diol 7-glucoside 
(S)-Naproxen 

(E)-Suberenol 


M/Z 


229.1222 
213.0916 
307.1186 
347.2211 


245.1168 


339.0543 
487.3413 
269.1745 
351.2157 


245.0122 


229.0866 
263.1276 
245.1170 
263.1274 
261.1118 
345.2421 
233.1170 
409.3462 
469.3305 
487.3424 
515.1190 
537.3193 
369.1917 
261.1131 
279.1228 
249.1117 
215.1066 
295.1534 
229.0866 


259.0974 


Retention 
time 
3.9183 
5.8296 
3.1032 
6.0994 


4.2623 


4.3228 
6.5588 
5.3919 
6.2181 


3.2813 


4.4459 
3.5650 
5.6671 
4.5786 
5.0895 
6.4991 
4.8241 
8.6280 
6.9737 
8.1855 
3.3125 
8.4688 
5.9223 
4.5816 
4.0935 
4.9490 
4.4806 
3.2376 
4.1531 


5.1608 


Formula 


CisH2004 
СНО 
CisH2006 
C21H3004 


CısH1603 


CısH1607S 
C30H46Os 
CisH24O4 
CooH3005 
CoHi1006S 


CisH16O4 
CisHisO4 
CisH16Os 
CisHisO4 
CisHi6O4 
CxH3203 
C14H16O3 
C29H440 
C30H4404 
C30H4805 
СНО 
C31H4806 
CisHosOs 
CisH200s 
CisHisOs 
CisH16O4 
СНО 
CicH2607 
C14H1403 


CisH16O4 


VIP_value 


6.1604 
4.4739 
4.4419 
3.8185 


5.543 


4.484 
6.0746 
5.4712 
4.1175 


4.4749 


4.328 
4.0153 
4.3166 
4.4066 
4.3101 
4.721 
4.5457 
4.1264 
3.7569 
3.7631 
3.6367 
3.7078 
3.6926 
3.8636 
4.0679 
3.9053 
3.9638 
3.9105 
4.7948 


4.7632 


P value 


1.659E-11 

1.438E-08 

6.684E-07 
0.000006894 


2.675E-10 


8.945E-08 
0.000002508 
2.113E-10 
0.001304 


2.69E-08 


6.039E-14 
0.000001038 
2.493E-13 
1.9E-10 
8.277E-11 
4.111E-12 
1.063E-09 
1.963E-07 
1.663E-08 
4.83E-09 
2.653E-08 
1.584E-09 
1.988E-16 
5.626E-08 
3.241E-12 
3.315E-11 
2.035E-08 
3.238E-08 
3.116E-16 


7.869Е-13 


